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Abstract. The Heine-Stieltjes Theorem describes the polynomial solutions, 
(v, /) such that T(f) = vf, to specific second order differential operators, 
T, with polynomial coefficients. We extend the theorem to concern all (non- 
degenerate) differential operators preserving the property of having only real 
zeros, thus solving a conjecture of B. Shapiro. The new methods developed are 
used to describe intricate interlacing relations between the zeros of different 
pairs of solutions. This extends recent results of Bourget, McMillen and Vargas 
for the Heun equation and answers their question on how to generalize their 
results to higher degrees. Many of the results are new even for the classical 
case. 



1. Introduction and Main Results 

Special cases of the following general eigen-problem have been studied frequently 
since the 1830's. 

Problem 1.1. Consider a linear operator T : C[z] — > C[z]. Describe all pairs 
(v, f) £ C[z] x C[z], with f 0, such that T(f) = vf. 

Recently the problem has received new interest from different perspectives, see 

e.g. m ei Gang. 

We say that v above is a (generalized) Van Vleck polynomial and that / is a 
(generalized) Stieltjes polynomial for T. The following theorem is known as the 
Heine-Stieltjes Theorem [15] and is due to Stieltjes [14], except for the statement 
regarding the location of the zeros of the Van Vleck polynomials in (3), which is 
due to Van Vleck [TB]. 

Theorem 1.2 (Heine-Stieltjes Theorem). Suppose that T : C[z] — > C[z] is a dif- 
ferential operator of the form T — Cj2(z)D 2 + Q\(z)D, where D — d/dz, Q%(z) — 

IT^iv^-ay) and 

d 

Qi{z) = Q2{z)T j ^^, 

* — ' z — a* 

j=l 3 

where a\ < ■ ■ ■ < ay are real, d > 2, and 71, . . . are positive. Let n > 2 be an 
integer. Then 

(1) There are exactly (™ + ^~ 2 ) pairs (v,f) G C[z] x C[z] such that T(f) = vf 
and f is monic of degree n; 

(2) For each Van Vleck polynomial there is a unique monic Stieltjes polynomial; 
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(3) The zeros of v and f in (1) are real, simple and belong to the interval 
(ai,ay). 

The problem studied in this paper is how to generalize the Heine-Stieltjes The- 
orem to larger families of linear operator. Which properties of T in the Heine 
Stieltjes Theorem are essential? We will see that the essential property of T is that 
it preserves the property of having only real zeros. 

A polynomial / £ M[z] is hyperbolic if all its zeros are real and a linear operator 
T : WL[z] — » R[z] preserves hyperbolicity if T(f) is hyperbolic or identically zero 
whenever / is hyperbolic. Such operators were recently characterized in [2]. Our 
extension of Theorem 11.21 is to differential operators, of finite but arbitrary order, 
that preserve hyperbolicity. Let 

N 

T=Y1 Qk(z)D\ where Qm(z)Q n (z) £ 0, (1) 

fc=Af 

be a differential operator with polynomial coefficients. The number 

r = max{deg Q k - k : M < k < N} 

is called the Fuchs index of T, and T is non- degenerate if degQjv = N + r. 

For a positive integer to, let [m] — {1,...,to}. An ordered partition of [to] 
into two parts is a pair (A, B) where A, B C [to], A U B = [to] and A fl B = 0. 
Let / and g be two hyperbolic polynomials of degree r and n, respectively and 
let ct\ < «2 < • • • < ctn+r be the zeros of fg. Let further (A, B) be an ordered 
partition of the set [n + r], such that \A\ = r and \B\ = n. The pair (/, g) is said 
to have zero pattern [A, B) if the zeros of / are {c^ : i £ A} and the zeros of g 
are {o^ : i £ B}. Note that the zero pattern is unique if and only if / and g are 
co-prime. 

It was conjectured by B. Shapiro 12. 13 , after extensive computer experiments, 
that a Heine Stieltjes Theorem should hold for non-degenerate hyperbolicity pre- 
serving differential operators. The solution to this conjecture, Theorem 11.31 is our 
main result. The reason why Theorem 11.21 is a special case of Theorem 11.31 is 
explained in Remark 12.31 

Theorem 1.3. Let T = J2k=o Qk(z)D k , where Q N {z) ^ 0, be a non-degenerate 
hyperbolicity preserving differential operator with nonnegative Fuchs index r, and 
suppose that n £ N is such that there is no 8 £ K such that Qq(0) = • • • = Q n {6) = 0. 
Then 

(1) If f\ifi are monic polynomial for which {v,fi) and (v,f2) are solutions to 
Problem \l.l[ then fx = fa; 

(2) For all SC. [n + r] of cardinality r there is a unique pair (v, f) of hyperbolic 
polynomials with zero pattern (S, [n + r] \ S) such that T(f) — vf, with f 
monic of degree n. These are all Stieltjes and Van Vleck polynomials. Hence 
there are exactly (™^ r ) pairs of polynomials (v, f) such thatT(f) = vf, with 
f monic of degree n; 

(3) v and f are co-prime and their zeros lie in the convex hull of the zeros of 
Q e , where e — mm(n, N); 

(4) If there is no € R such that Qi(8) = ■ ■ ■ = Q n {0) = then v and f have 
simple zeros which all lie in the interior of the convex hull of the zeros of 
Q e , where e = min(n,iV). 
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That v and / are co-prime in the original Heine-Stieltjes theorem is due to Shah 

It follows from Lemma [2.71 below that the coefficients Qj(z), for < j < N, 
in Theorem 11.31 are hyperbolic. Hence the convex hull of the zeros of Q e (z) in 
Theorem II .31 (3) and (4) is a closed interval. 

Example 1.4. A natural class of differential operators satisfying the hypothesis of 
Theorem 11.31 consists of operators of the form 

r (/)0) = D N ~ M {PD M f){z) = ( N k lM) PiN ~ k) ^ Dkf ^ 

k=M ^ ' 

where P is a hyperbolic polynomial of degree N + r, with no zero of multiplicity 
larger than N — M. 

The following theorem offers a way of (in theory) generating Stieltjes and Van 
Vleck polynomials. Note that there is no restriction on the common zeros of the 

Theorem 1.5. Let T = £f =0 Q k (z)D k be a non- degenerate hyperbolicity preserv- 
ing differential operator with nonnegative Fuchs index r. Suppose that n £ N is 
such that T(z n ) ^ and that A C [n + r] is a set of cardinality r. Let a be the 
smallest zero of Q e {z), where e = min(n, N). 

Define two sequences of polynomial {.fi}°^ and {vi}°°^ 1 recursively by fo{z) = 
(z — a) n and 

T{fi) = Vi+ifi+i, 

where Wi+i/^+i is the unique factorization ofT(fi) such that fi + \ is monic and the 
pair (ui-fij fi+i) has zero pattern (A, B), where B = [n + r] \ A. Then 

lim fi = f and lim Vi = V, 

i — >oo i — too 

where T(f) = vf and (v,f) have zero pattern (A,B). 

In Section 2] we use Theorem 11.51 to derive interlacing properties of the zeros of 
different Stieltjes and Van Vleck polynomials. These generalize recent results of 
Bourget, McMillen and Vargas [5] and Bourget and McMillen ;4] derived for the 
Heun equations, and answers the question raised in [4] on how to generalize their 
results to Van Vleck polynomials of higher degree. 

2. Properties of Hyperbolicity Preservers 

Here we collect the results on hyperbolicity preserving differential operators 
needed to prove Theorem 11.31 Some results are new and of independent interest. 
A polynomial f(z\, . . . , z n ) £ C[zi, . . . , z n ) is stable if f{z\, . . . , z n ) ^ whenever 
Im(zj) > for all 1 < j < n. Hence a univariate polynomial with real coefficients is 
stable if and only it is hyperbolic. The following characterization of hyperbolicity 
preservers in the Weyl algebra was given in pQ . 

Theorem 2.1. Let T = J2k=M Qk{z)D k , where Q k {z) e R[z] for all M <k < N 

and Qm(z)Qn{z) ^ 0. Let further Gt{z, w) — Ylk=M Qk(z)w N ~ h £ M.[z, w] . The 
following are equivalent: 

(1) T preserves hyperbolicity; 

(2) Gt{z,w) is stable; 
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(3) Gt(z,w) can be expressed as 

(z, w) = ± det(zA + wB + C) 

where A and B are positive semi-definite symmetric matrices and C is 
symmetric. 



Remark 2.2. Using Theorem 12.11 one can easily construct linear operators T for 
which Theorem ll.3l is applicable: Choose A, B, C to be symmetric (N + r) x (N + r) 
matrices satisfying: A is positive definite, B is positive semi-definite and has as 
an eigen- value of multiplicity r, and C is such that det(0A + Bw + C) ^ for all 
9 e K. Then choose T such that G T (z, w) = ± dct(zA + wB + C). 

Remark 2.3. Let T be as in Theorem 11.21 Suppose that Gx[z,w) — Q-2(z) + 
Qi(z)w — where Im(z) > 0. Since Qi is hyperbolic and 



we have w ^ and 



L,,,^!.. Mf> <0 

\z-ajY 



Im( W - 1 )=5: 7 ,^%>0. 
z — ' z — a, r 



Hence Im(w) < 0. Thus Gt(z,w) is stable and by Theorem 12.11 T preserves 
hyperbolicity. Note also that Qi and Qi are co-prime. Therefore Theorem 11.21 
follows from Theorem 11.31 In fact, the linear operators considered in Theorem 1 1.21 
are precisely those with M = 1 and N — 2 in Theorem 11.31 

Proposition 2.4. Let T = J2k=M Qk{z)D k , where Q k {z) £ R[z] for all M < k < 
N and Qm{z)Qn{z) ^0, be a hyperbolicity preserver. Let r = degQjvf — M. Then 

(1) dcg Q t < r + i for all M < i < N; 

(2) If ai is the coefficient in front of z r+l in Qi(z), then all the zeros of the 
polynomial Y]^-m a-iZ 1 are real and nonpositive. In particular, the nonzero 
aiS have the same sign, and the indices i such that a.i =/= form an interval; 

(3) T(z n ) is of degree n + r for all n > M and the leading coefficient of T(z n ) 
is equal to \ n = i!ai(™). In particular, if am > 0, then 

An+i — A 



i 

for all n>M, unless M = N = 0. 

Proof. Property (1) is a direct consequence of Theorem 12.11 and the fact that the 
support of a stable polynomial forms a jump system, see [51 Theorem 3.2]. 

Let 7 > 0. Then -y- r - N G T (jz, 7 u>) is stable by Theorem O By Hurwitz' 
thcorenu so is the homogeneous polynomial 

N 

g(z,w) = lim 7~'~ Gri'JZjjw) = > aiz r+t w ~ l . 

7 — >oo ^ — 4 

i=M 

A property of homogeneous stable polynomials is that all nonzero coefficients have 
the same phase, see Theorem 6.1]. The polynomial g(z, 1 + ie) is stable if e > 0. 
Hence, by Hurwitz' theorem, g(z, 1) is stable and hence hyperbolic. Since the 



^See [7] p. 96] for an appropriate multivariate version. 
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nonzero coefficients of g(z, 1) have the same sign, the zeros of g(z, 1) are nonpositive. 
The indices i such that ai ^ of such a polynomial form an interval. 

Property (3) follows easily from (1) and (2). □ 

Let ai < OL2 < • • ■ < a n and (3\ < 02 < ■■■ < 0m be the zeros (counted 
with multiplicities) of two hyperbolic polynomials / and g, with deg/ = n and 
degg = to. We say that these zeros interlace if they can be ordered so that either 
Oil < Pi < ct2 < /?2 < • • • or (3\ < ai < /?2 < a 2 < • • • j in which case |n — m| < 1. 
Note that by our convention, the zeros of any two polynomials of degree or 1 
interlace. It is not difficult to show that if the zeros of / and g interlace then the 
Wronskian W[f, g] ■= f g — fg is either nonnegative or nonpositive on the whole 
real axis E, see, e.g., [TU]. In the case that W[f, g] < we say that / and g are in 
proper position, denoted / <C g. For technical reasons we also say that the zeros 
of the polynomial interlace the zeros of any (nonzero) hyperbolic polynomial and 
write <C / and / <C 0. Note that if / and g are (nonzero) hyperbolic polynomials 
such that / <C g and g <C / then / and g must be constant multiples of each other, 
that is, W[f,g] = 0. 

The following theorem is a version of the classical Hermite-Biehler theorem, see 

Theorem 2.5 (Hcrmitc-Bichlcr theorem). Let h := f+ig £ C[z], where f,g £ M.[z]. 
Then h is stable if and only if g <C / ■ 

The next proposition follows from the fact that a hyperbolicity preserving dif- 
ferential operator also preserves stability when acting on complex polynomials, see 

P 

Proposition 2.6. Let T be a hyperbolicity preserving differential operator. Suppose 
that f,g 6 M.[z] are such that f -C g. Then T(f) <C T(g), unless T(f) — T(g) = 0. 

Lemma 2.7. Let T = J2k=M Qk(z)D k , where Q k {z) G R[z] for all M < k < N 
and Qm{z)Qn(z) ^ 0, be a hyperbolicity preserver. Then Qj <C Qj+i for all 
M < j < N—l. Ln particular Qj is hyperbolic or identically zero for all M < j < N . 

Proof. Let S be the linear operator on M.[z, w] defined by S(z k w i ) = z k w e if 
£ € {N - j,N - j + 1} and S{z k w e ) = otherwise. It follows fr om the theo- 
rem characterizing multidimensional multiplier sequences in [1] that S preserves 
stability. The polynomials Gt(z,w) in Theorem 12. II is stable and 

S(G T ) =w N - j (Q j (z) + wQ j -i(z)). 

Hence Qj(z) + wQj-i(z) is stable and so is the polynomial Qj(z) + iQj-i(z). The 
lemma now follows from the Hermite-Biehler theorem. □ 

If / is a hyperbolic polynomial we let /(/) denote the smallest interval containing 
all the zeros of /. 

Lemma 2.8. Let T = Ef=M Qk(z)D k be a non- degenerate hyperbolicity preserver, 
where QmQn ^ 0. Suppose that f is a hyperbolic polynomial of degree n > M. 
Then 

I(T{f)) C I(Q e f), 

where e = min(n, N). 
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Proof. Since T is non-degenerate we know by Proposition ^. 4l that deg Qk = deg Qm + 
k, for all M < k < N. By Proposition 12.41 the sign of the leading coefficients of 
the Qfe's is the same for all M < k < N . Since Qm -C Qm+i *C • • • <C Qn, by 
Lemma [271 it follows that I{Qif {i) ) C 7(Q e /) for all M < i < n. Hence, all the 
polynomials Qif^ are nonzero and have the same sign outside I(Q e f). Since T(/) 
is a nonnegative sum of these polynomials we have 7(T(/)) C I{Q e f). □ 



The next lemma is often referred to as Obreschkoff 's theorem, for a proof see 
[TO] or [3J. 

Lemma 2.9. Lei /,g 6 R[z]. TTien aZZ nonzero polynomials in the space 

{af + f3g : a, /3 e K} 
are hyperbolic if and only if f <C 5, 5 <C / or f = g = 0. 

Lemma 2.10. Lei F be a hyperbolic polynomial of degree d and let f be a hyperbolic 
polynomial. Suppose that (z - 6) N | F{D)f, where N > 2. Then (z - 6) N+d \ f or 
F{D)f = 0. 

Proof. The lemma follows from the case d — 1 by induction on d. Hence let 
F = a + z where a G M. It follows from e.g. Lemma l2~9l that the zeros of / and the 
zeros of /' interlace the zeros of af + f . Hence (z — 6') 7V_1 | /' and (z — 8) N ~ 1 \ f, 
and since N > 2 we have (z - 9) N \ f. Since /' = F(D)f - af, (z - 6) N \ F(D)f 
and (z -9) N | / we have that (z - 9) N \ f . It follows that (z - 0) N+1 \ f as was 
to be shown. □ 



Lemma T2.10I generalizes to differential operators as follows. 

Proposition 2.11. LetT ~ Ylk=M Qk(z)D h , where QmQn ^ 0, be a hyperbolicity 
preserver and let f be a hyperbolic polynomial with n = deg / > M. Suppose that 
(z-9) d | T{f), where d>2. Then either Q M {9) = Q M+ i{9) =■■■ = Q n (9) = or 

(z-9) d + t \f(z), 



where t = max{fc < N : Q k {9) ± 0}. 



Proof. Let / and 9 be as in the statement of the proposition and let 

A " T(/( fc ))(0) 



G(8,w)= J2 Qk{0)f {k \w + 9) = Y J 



fc! 

k=M k=0 



-w k . 



If Qm(Q) = Qm+i{9) = ... = Q n {9) = there is nothing to prove so assume 
otherwise. It follows from Theorem 12.11 and the fact that stability is preserved 
under setting variables to real values (see e.g. [3]) that the polynomial F(z) = 
ELm Qk(0)z k is hyperbolic. By hypothesis F(D)[f(w+9)} =G(9,w) ^0. Clearly 
the zeros of f^> and f( l+1 > interlace for all < i < n. By Lemma [2.61 the zeros 
of and T(/( i+1 )) interlace for all < i < n - M. Since (z - 9) d \ T(f) it 

follows that - 9) d - 1 | T(f^) for all < i < d. Hence w d \ G(9,w) which by 
Lemma [2T0] implies w d+t \ f{w + 6), that is (z - 9) d+t \ f(z). □ 
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3. Proofs of the Main Results 
We will make use of the following generalization [12] of Heine's theorem [8] . 

Theorem 3.1 (B. Shapiro). Let T be a non- degenerate differential operator as in 
([T]). and let neN. Suppose that T has Fuchs index r > 0. Then there are at most 
( n ^ r ) distinct Van Vleck polynomials v, for which T(f ) — vf for some f <E C[z] of 
degree n. 

Lemma 3.2. Let T = J2k=M Qk(z)D k , where Q M Q N ^ and M < N, be a 
hyperbolicity preserver with Fuchs index r > 0. Then there is exactly one monic 
Stieltjes polynomial to each non-zero Van Vleck polynomial. 

Proof. The set of Stieltjes polynomials (union the zero polynomial) corresponding 
to the same Van Vleck polynomial, v, is a linear space. Hence if there were two 
monic Stieltjes polynomials corresponding to v, then there would be two Stieltjes 
polynomials of different degrees m > n, also corresponding to v. However, by 
Proposition 12. 4i the leading coefficient of a Van Vleck polynomial pairing with a 
Stieltjes polynomial of degree k is and since {\i}°°^ M is strictly increasing we 
have a contradiction. □ 

Proof of Theorem \1.3{ Let T and n be as in the statement of the theorem and let 

5 C [n + r] be of cardinality r. We first prove that there is a pair (v, f) of hyperbolic 
polynomials with zero pattern (S, [n + r]\S) such that T(f) = vf with / of degree 
n. Let I be the convex hull of the zeros of Q e . Hence 7 is a closed bounded 
interval, possibly without interior. Define a function <f> : I n — > I n as follows. If 
a = (pti, . . . , a n ) 6 7™ form the polynomial h(z) — " a j)- By Proposition 
12. 4\ T(h) is hyperbolic of degree n + r. Hence we may write T(h) uniquely as vg 
where g is monic and (v, g) has zero pattern (S, [n + r]\S). Let 4>(a) = (J3\, . . . , /3„), 
where /3± < • • • < /3 n are the zeros of g. By Hurwitz' theorem <fi is continuous, which 
by Brouwer's fixed point theorem implies that <fi(a) = a for some a £ I n . Hence 
T(f) = vf for some (v, f) with zero pattern (5, [n + r] \ S). 

Next we prove that if (v, f) is a pair of hyperbolic polynomials with zero pattern 
(S, [n + r] \ S) such that T(f) — vf with / of degree n, then v and / are co-prime. 
Assume otherwise that f(0) = v(9) — and that the multiplicity of 9 in / is m. 
Then (z - 6) m+1 | T(f) which by Proposition EU implies (z - 6) m+1 \ f, which is 
a contradiction. Note that this means that the zero pattern of (v, f) is unique. 

(1) is just a special case of Lemma T3. 2 1 

From what we have proved so far we deduce that there are at least ( n ^ r ) distinct 
Van Vleck polynomials. By Theorem 13.11 these are all Van Vleck polynomials, and 
by the above there is a unique Stieltjes polynomial to each Van Vleck polynomial. 

It remains to prove (4) for a given pair (v, f) of Van Vleck and Stieltjes polyno- 
mials. Suppose that / has a zero of multiplicity m > 1. Then [z — 9) m \ T(f). 
By Proposition 12.111 (z — 9) m+1 | /, which is a contradiction. If v has has a zero 

6 of multiplicity m > 1, then (z - 6) m \ T{f). Hence (z - 6) m+1 \ f by Proposi- 
tion [nHU This contradicts the simplicity of the zeros of / just derived. Let a be 
an endpoint of /. Since all zeros of / are simple and located in I we know that 
f'(a) ■ ■ ■ f( n \a) 7^ 0. As in the proof of Lemma [2.81 we note that the signs of all 
nonzero Qi(a)f^(a), for < i < e, are the same. By hypothesis Qi(a)f^(a) ^ 
for some i, which proves that T(/)(a) ^0. □ 
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Proof of Theorem I J. 51 Let I and (j> : I n — » I n be as in the proof of Theorem 11.31 
Let further < denote the usual product partial order on R™. We claim that if 
a, (3 G K n are weakly increasing, i.e., a± < ■ ■ ■ < a„, (3i < • • • < /3„ and a < j3, 
then 0(a) < 0(/3). To see this write a <C /3 if ai < /?i < 0:2 < $2 < • • • < ct n < (3 n . 
Clearly a < [3 whenever a -C [3. By Proposition 12.61 it follows that (j)(a) <C (f)(0) 
whenever a, (3 S K™ are weakly increasing and a <C However, if a < /3 are 
weakly increasing then there is a sequence, {7 fe }fcLo ^ sucn that a = 7 <C 
7 1 < • • • < 7 m = /3. Hence 0(a) < 0(7 1 ) <C • • • < 0(/3), and thus (f>(a) < 4>((3) 
which proves the claim. 

Define a sequence {a fc }^L C I n recursively by a = (a,..., a), and a k+1 = 
<p(a k ) for all k G N. Thus the coordinates of a k are the zeros of Since a is 
the unique smallest element in T" with respect to <, and 4> : I n — > 7™, we have 
a < a 1 . Since preserves the partial order < we have a < a 2 < a 3 < • • • . By 
compactness, the sequence {a k }^L converges to a limit a = (aj., . . . , a„) £ I n . It 
follows that a is a fixed point of <fi and that lirm^oo fi = f := JlILiC 2 — Since 
all the ViS have the same leading coefficient, A„, it follows that lim^oo Vi = v and 
that T(f) = vf where (v,f) has zero-pattern (A, B). 

□ 

4. Interlacing Properties 

In this section we will prove interlacing properties of zeros of Stieltjes polynomial 
and of Van Vleck polynomials. The first two theorems concern Stieltjes polynomi- 
als of the same or consecutive degrees. These theorems generalize the theorems 
obtained by Bourget and McMillen [4] for the Heun equation. In [4] it was asked 
if and how their results could be extended to the case of Fuchs index greater than 
one. Theorems 14.31 and 14.51 answers this. Our theorems concerning Van Vleck 
polynomials generalize to arbitrary nonnegative Fuchs index the results of Bourget, 
McMillen and Vargas [5] regarding the Heun equation. 

Definition 4.1. Let ( Jl ^) denote the set of all n element subsets of [n + r], and 

let T — J2k=M Qk(z)D k , where QmQn ^ 0, be a non-degenerate hyperbolicity 
preserver with nonnegative Fuchs index r. Let further B E (^ r ') and A — [n + 
r] \ B, where n > M. We denote by (va,ti, /s,n) the pair of polynomials with 
zero-pattern (A,B) afforded by Theorem 11.51 That is, fs.n is monk of degree n 
and T(f B , n ) = v A , n fB,n- 

Henceforth it will tacitly be assumed that n > M as in Definition 14. II 
A pair (x, y) of real vectors x = (xi, . . . , xt) and y = (y 1: . . . , y^) are said to be 
in proper position, written xCy, if 

x\ < yi < x 2 < V2 < ■ ■ ■ < x k < y k . 

If A is a set of integers let [A] denote the vector obtained by ordering its elements 
increasingly. Define a relation — > on (^ r ') by: A — > B if 2[A] < 2[B] + 1, where 
1 = (1,...,1). 

The relation — > may be described differently as follows. If x = (x\, . . . , £„+,.) G 
M. n+r and B G ('™^ r '), let be the vector obtained by deleting the coordinates 
labeled by [n + r]\B. Then B -> C if for all x, y G W l+r with x < y we have 
xb < yc- 
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Lemma 4.2. Suppose that the polynomials f and g have interlacing zeros and 
that 9 G K is a common zero of f and g, of multiplicity one in each. Then 9 is 
a zero of multiplicity exactly two of the polynomial g'{9)f{z) — f'(9)g(z), unless 
g'(6)f(z)-f(6)g(z) = 0. 

Proof. Let F{z) = g'{9)f(z) - f'{9)g{z) and suppose that F(z) ^ 0. Clearly 
F(9) = F'(9) = 0. By Lemma 12.91 the zeros of F(z) and f(z) are interlacing. 
Hence 9 is a zero of F (z) of multiplicity at most two. □ 

The following theorem describes interlacing relationships betwen Stieltjes poly- 
nomials of the same degree. 

Theorem 4.3. Let T be a non-degenerate hyperbolicity preserver with nonnegative 
Fuchs index r and let B,C G ('"^ r ') ■ If B — > C then fs.n "C fc,n- 

Moreover, if there is no 9 G K for which Q\{9) = ■ ■ ■ = Q n (9) = 0, then fB,n 
and fc,n are co-prime, i.e., their zeros strictly interlace. 

Proof. Let {f l B n }?Z and {fc n }i^o be the sequences of polynomials constructed 
as in Theorem 11.51 for fB, n and fc,n, respectively. We prove by induction that 
fh.n ^ fc n f° r au integers i > 0. The induction start is trivial since f% n = fc n- 

Suppose that f l Bn fhn- Then, by Proposition 12.61 we have T(f Bn ) <C 
T(/p n ). The polynomials f B + n and f%t n are obtained by keeping the zeros of 
T(fhn) and T(fbn) that are labelled by B and C, respectively. Hence, by the 
(alternative) definition of — > we have f B + n <C fj^n' 

The first part of the theorem now follows by letting i — > oo. 

To prove the second part assume that fB,n{j) — fc,n{l) = for some 7 G M. 
By Theorem 1 1.31 fB,n and vb 1 ,n are simple-rooted and co-prime. The same holds 
for fc,n and vc,n- By Lemma 14. 2 [ 7 is a zero of multiplicity exactly two of the 
polynomial T(f B , n ~ afc,n) Sb fc,n, where 

a _ {vB',nfB,n)'in) 
(vC\nfc,n)'h) ' 

But then, by Proposition 12. 1 H 7 is a zero of multiplicity at least three of the 
polynomial fs,n — afc,n-, which contradicts Lemma 14.21 □ 

Remark 4.4. Given B = {b\ < ••• < b„} C [n + r] there are exactly 61(62 — 
61) • • • (b n — b n -i) sets A G ('"^ r ') for which A — » B. Indeed, to construct A we 
pick an integer from the interval [6^ + l,6fc+i], independently for each k. By a 
counting argument there are exactly ( 2n+r ) relations (including the trivial ones) in 

cr 1 )- 

Similarly, a pair (x, y) of real vectors x — (xi, . . . , Xf.) and y = (yi, . . . , yu+i) 
are said to be in proper position, written x -C y, if 

Vi < xi < 2/2 < x 2 < ■ ■ ■ < x k < yk+i- 

The relation -> extends naturally to concern ( [n + r] ) x ( [ "^+{' rl ): A ->■ B if 2[A] <C 
2[B} + 1. Again, A —> B \i and only if for all x G K Il+r and y G U n + r + 1 with x < y 
we have x^ <C ys- 

The next theorem describes interlacing relationships between Stieltjes polyno- 
mials of consecutive degrees. Its proof is very similar to that of Theorem 14.31 and 
is therefore omitted. 
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Theorem 4.5. Let T be a non-degenerate hyperbolicity preserver with nonnegative 
Fuchs index r and let B G ( [n + r] ) and C G ( [n ++j r] ) • If B — > C then f B , n < fc,n+l- 
Moreover, if there is no 9 S R for which Q\(9) = ■ ■ ■ = Q n (9) = 0, then fs,n 
and fc.n are co-prime, i.e., their zeros strictly interlace. 

We have the following theorem describing interlacing properties of Van Vleck 
polynomials that have Stieltjes polynomials of the same degree. 

Theorem 4.6. Let T be a non-degenerate hyperbolicity preserver with nonnegative 
Fuchs index r and let A G ([" + ' r ~ 1 ]) . Then va n *C va+i n , where A+l := {a+1, a G 
A}- 

Moreover, if there is no 9 G R for which Q\(9) = ■ ■ ■ = Q n (9) = 0, then VA,n 
and VA+\.n are co-prime, i.e., their zeros strictly interlace. 

Proof. Clearly (A+l)' -> A', where (A + l)' = [n + r]\(A + l) and A' = [n + r]\A. 
By Theorem 14.31 f(A+i)' *C /a' and by Proposition ^. 61 

VA+if(A+iy = T(/( A+1 )/) < = iu/a'- 

Let Xi < • • • < x n + r and yi < ■ ■ ■ < y n +r be the zeros of VA+\f(A+\)> an d VAfA>, 
respectively. Then, for nei, 

a:i < 2/i < ■ ■ • < x a < Va < x a +i < y a +i <'■■< x n+r < y n+r . 

Since the zeros of va are {y a ■ a G ^4} and the zeros of va+x are {x a+ i : a G A}, it 
follows that va *C va+i- 

The proof of the strict interlacing of the zeros follows just as in the proof of 
Theorem E2 □ 

The next theorem describes interlacing properties of Van Vleck polynomials that 
have Stieltjes polynomials of consecutive degrees. The proof is very similar to that 
of Theorem 14.61 and is therefore omitted. 

Theorem 4.7. Let T be a non-degenerate hyperbolicity preserver with nonnegative 
Fuchs index r and let A G ('™^ r '). Then va.u+i *C va.u *C t>4.+i,n+i- 

Moreover, if there is no 9 G R for which Q\(9) = ■ ■ ■ = Q n (9) = 0, then the 
interlacing of the zeros is strict. 

Suppose that T — Y1%=m Qk(z)D k is a non-degenerate hyperbolicity preserver 
with Fuchs index r = 1, and that Qm(9) = ■■■ = Q n (9) for no 9 G R. Then 
the Van Vleck polynomials are of degree one and there are exactly n + 1 of them. 
Define the nth spectral polynomial, p n (z), to be the monic polynomial whose zeros 
are precisely the zeros of the n + 1 Van Vleck polynomials. 

Corollary 4.8. Suppose thatT = Qk( z )D k is a non- degenerate hyperbolic- 

ity preserver with Fuchs index r = 1, and that Qm(9) = ■ ■ ■ = Q n (9) for no 9 G R. 
Then the zeros of p n and p n +i interlace. Moreover p n and p n +i have no zeros in 
common. 

Proof. The interlacing property follows immediately from Theorem 14.71 That p n 
and Pn+i have no zeros in common is a consequence of Lemma 13.21 □ 

For the Heun equation Corollary [48] reduces to the main result, Theorem 2.1, of 

0. 
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5. Open Questions and Further Directions 

The Jacobi orthogonal polynomials arise as the Stieltjes polynomials correspond- 
ing to the case when Q^iz) is of degree 2 in the classical Heine-Stieltjes theorem. 
Given a hyperbolicity preserving differential operator as in Theorem 11.31 (4) of 
Fuchs index r — 0, it is natural to ask whether the Stieltjes polynomials form an 
orthogonal family. This was suggested by T. McMillcn. 

The linear operators considered in this paper are all finite order differential op- 
erators. One may try to extend our results to hyperbolicity preservers that are not 
finite order differential operators. There is a natural class of such operators for 
which there are pairs of Van Vleck and Stieltjes polynomials for each possible zero 
pattern. The problem is to ensure the uniqueness of the Stieltjes- Van Vleck pair 
given the zero pattern. 

Acknowledgments. I thank Julius Borcea and Boris Shapiro for introducing me 
to the Heine-Stieltjes problem, and Julius Borcea for bringing [4] to my attention. 
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